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ABSTRACT 

We have studied the relation between the structural relaxation and the frequency dependent 
thermal response or the specific heat, c p (u), in a supercooled liquid. The Mode Coupling 
Theory (MCT) results are used to obtain c p (u) corresponding to different wavevectors. Due 
to the two-step relaxation process present in the MCT, an extra peak, in addition to the low 
frequency peak, is predicted in specific heat at high frequency. 



Introduction 



Understanding the complex relaxation behavior in supercooled liquids has been a field of 
much research interest in recent times. In this regard the response of a system to an energy 
fluctuation namely the frequency dependence of specific heat of a supercooled liquid has been 
investigated by a number of authors. Generally, specific heat is a property that is usually 
linked to the thermodynamic property of a system . The pioneering experiments done by 
Birge and Nagel [0 have studied the dynamic response in glassy systems namely glycerol 
and propylene and obtained interesting dynamical response behavior, expressed in terms of 
a frequency dependent specific heat 0. In an experiment usually the frequency dependent 
product of thermal conductivity(K) and specific heat, kc p {uj) is measured. However, in the 
temperature range 190°K — 220°fT, over which we are interested here, k has weak frequency 
dependence 0,0 and hence the dynamics observed in the product KCp{u) is entirely due to 
the frequency dependence of the specific heat. The theoretical modeling for the frequency 
dependence of specific heat in a supercooled liquid has been studied by various authors 
0,0. New internal mode was proposed to be present in the supercooled liquid to explain 
the frequency dependent response. In a simple analysis, Zwangig however had argued that 
the frequency dependence of the specific heat can be obtained without introducing any such 
internal mode. This work showed that what is measured as the frequency dependent specific 
heat is actually related to that of the longitudinal viscosity in the liquid. In this model the 
dynamics of fluctuations around the equilibrium was studied in terms of a simple set of slow 
variables of Hydrodynamics of fluids. These equations of motion used were the conservation 
laws of mass, momentum and energy in the system. The resulting formula for the specific 
heat is equivalent to linking of structural relaxation in a supercooled liquid to the frequency 
dependence of the specific heat. In the present work we take the data from the Specific heat 
measurement |l[ and extract the frequency dependence of the viscosity as will be required 
from such a formulation proposed in Ref. 0. We then address the question if this value 
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of the longitudinal viscosity will indeed be self consistent with independent measurements 
on the structural properties. The basic idea is to consider the frequency dependence of the 
specific heat solely in terms of the structural relaxation in the supercooled liquid. 

For structural relaxation - the microscopic model for the liquid dynamics namely the 
Mode Coupling Theory (MCT) has been studied by a number of authors in recent years 
[§. § . A simple application of the MCT to fit the specific heat data indicated a very large 
exponent is required to match the experimental data to power law divergence. In the present 
work we also use the MCT as a model for structural relaxation and obtain the corresponding 
frequency dependence of the longitudinal viscosity. We then use it to predict the behavior 
of the specific heat with frequency as predicted from the theory proposed in Ref 0. In 
the microscopic model of the Mode coupling Theory the wave number dependence of the 
longitudinal viscosity is obtained using proper input for the structure factor of the liquid. For 
this purpose standard results from the integral equations for simple liquids are used for the 
structure factor. The frequency dependent specific heat is then computed for different wave 
numbers. Thus the effect of the response to heat fluctuations can be computed over different 
length and time scales in the present approach. While this extends the theory with scope of 
further comparison, the main goal of the present work is to test if the frequency dependence 
of the specific heat can be understood solely in terms of the structural relaxation and if 
the two sets of measurements agree in a self-consistent manner. The paper is organized as 
follows. In the next section, Sec. II, we consider the schematic model for the time dependence 
of the viscosity and, in Sec. Ill, we compare the theoretical results with the experimental 
observations. A wavenumber dependent calculation for the specific heat is presented in the 
Sec. VI. In Sec. V, we present the Mode Coupling results for the specific heat over different 
length-scales and temperatures. In the last section we discuss the results. 
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II. Frequency Dependent Specific Heat 

Since we are concerned here with the dynamic properties of a supercooled liquid, an obvious 
choice is to consider a hydrodynamic model for the system. To start with, we write down 
the linearized hydrodynamic equations, 

p P + V.g = (1) 

?2L + ViP _ ^ - {\v + C)V(V.&) = (2) 

^4 5T+KV2T+ ^(f)/-^°' (3) 

governing the time evolution of fluctuations of conserved variables- mass density (p), momen- 
tum density(g) and the temperature(energy)T. Here p Q and T represent equilibrium density 
and temperature respectively and 5p and 5T are the fluctuations from the equilibrium val- 
ues. c v is the specific heat per unit mass at constant volume, 7] and ( are the shear and 
bulk viscosities respectively. The viscosity coefficients here are divided by the density. The 
fluctuation of the pressure P around the equilibrium value can be expanded to the lowest 
order in density and temperature as, 

r 

where we have assumed that the change of the Pressure functional with the density function 
at the equilibrium can be replaced by the equilibrium partial derivative - replaced by the 
corresponding thermodynamic quantity. Using the above equations the energy conservation 
equation @ reduces to the Fourier heat law for thermal fluctuations, 

luj5T = p(oo)V 2 6T, (5) 

with the frequency dependent thermal diffusivity p(oo) defined in terms of the specific heat 



The specific heat c p (uj) is expressed in the form, 

Kt(0) 

c p {u) =c v + {c p -c v ) . (7) 

A^l £J ) 

The quantity Kt{oj) is called the generalized bulk modulus and is given by 

Kt(uj) 

-^l = l + uvr(u), (8) 

and is expressed in terms of the reduced form T(u) = t]i{uj)/c? of the frequency dependent 
longitudinal viscosity rji(uj). In equation (§), K is the uj = limit of Kt(uj). Obviously 
for the liquid state with finite viscosity the zero frequency limit of Kt(oS) relates to the 
thermodynamic property of the supercooled liquid. The sound speed c Q is given by 

- f £ ) = t • (9) 

\ d P/ T Po 

A frequency dependent longitudinal modulus M(u), the inverse of compliance, is defined 
along a similar line as, 

^ = 7 +*aV), do) 

where 7 is the ratio of the long time limit of the specific heat c p {uj = 0) to c v . Equation (^) is 
the key formula used in this paper for testing the idea of modeling the frequency dependence 
in the specific heat solely in terms of the structural relaxation. In obtaining equation ([5]) 
one also needs to assume that the following self-consistent relation holds, 

Here we have expressed M in the dimensionless form as, M{uj) = M(u)/K . v = Cg/ fi D , 
with fi D = n/(p c v ) is the bare thermal diffusivity. We test the validity of the assumption 
([TT]) in the frequency range where the analysis with respect to experimental data is made. 
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III. Comparison of Experimental Data 

In this section we test the self-consistency in expressing the frequency dependent data on 
specific heat and Structural relaxation. For the supercooled liquids the relaxation over 
longest time scales, i.e. the a-relaxations, follows the stretched exponential behavior, 



rj(t) = rj exp 



T 



(12) 



where rj is the amplitude and j3 is the stretching parameter which defines the degree of de- 
viation from the exponential decay. In fitting the specific heat data we use the dimensionless 
form for the specific heat ratio, 

(7-1) 



c p (uj) — c v 



1 + 



(13) 



1 + io?r(a?) 

which reduces the formula in a dimensionless form. We fit the specific heat data of Ref. [TJ to 
the formula ( |T3"D using a simple stretched exponential ( |T2| ) relaxation function. In Fig. 1(a) 
and 1(b) we show the respective fitting of the experimental data for the real and imaginary 
parts of Cp(u) in the supercooled glycerol for three different temperatures, T = 201. 4°K, 
203.9°fC, 211 A°K. The arrows in Fig. 1(b) indicate the peak positions in the imaginary 
parts of the viscosities at the corresponding temperatures. In calculating the specific heat, 
the three parameters T(t = 0)), the relaxation time r and the stretching exponent (3 are used 
as the free parameters. 7=1.86. Using the best fit values of the parameters with the specific 
heat data we compute the structural properties of the liquid given by Modulus M(uj) defined 
in eqn. (ffH) and the longitudinal viscosity. The resulting behavior for these quantities are 
compared with the experimental results as shown in Figures 2 and 3. 

In Figure 2, we show the viscosity r\ in the zero frequency limit in units of K r where r D 



is the unit of time used. In the Inset we show the corresponding experimental data [|K| [11 
for the viscosity. The experimental data shown here is over a much wider temperature range 
(317°K ~ 190°i^) - both theoretical and experimental data agree with the Vogel Fulcher 



fit ( r] ~ exp[A/(T - T )] for T Q = 128°K and A = 2480°K shown in both the figure and 
the inset as solid lines. The zero frequency modulus K Q is roughly temperature independent 



1^1 over the range considered here. We notice that the viscosity increases by four orders 



of magnitude as the temperature is decreased over a small range (200°/^ — 220° K) near the 
glass transition temperature(T g =190° K). We define a normalized longitudinal modulus 

^ = MM-M(O) 

V ; M(oo)-M(0) V ; 

In Fig. 3 (a) and (b) we show the real and imaginary parts of M respectively denoted 
by M'(uj) and imaginary M"(u) against the frequency for the three different temperatures, 
T=203.9 o ir(solid), 211A°K (dashed), 221.5°K (dotted) lines. The frequency in each case is 
expressed in terms of the ratio with the corresponding peak position (u p ) in the imaginary 



part. The corresponding results from measurements on the modulus M [13] are also shown 
with filled circle. In Fig. 4 we show the plot of the peak positions as found in the fitting 
with different temperatures. As the temperature is decreased, peak in the imaginary part of 
the specific heat shifts towards the lower temperatures, signifying the slower relaxations in 
the system. The solid line indicate V-F fit with T Q = 128°K. In Fig 5 we show the frequency 
dependent Specific heat and the viscosity function ( in the inset ) at the same temperature. 
The peaks appear nearly at the same position on the frequency scale for the two quantities 
showing that the dominant time scales are same in the two cases. Finally we test the validity 
of the assumption (|ll|) that is crucial in reaching the Fourier heat law ([5]) with the frequency 
dependent specific heat - defined above. For this we calculate both the real(A'(co>)) and 
imaginary (A"(uj)) parts of (A(uj)) for the supercooled glycerol. In Fig. 6 we plot both the 
real and imaginary parts of A (a;) on the frequency range over which specific heat (frequency 
dependent) is observed. These figures clearly show that the quantity A(u) is much smaller 
as compared to unity over this frequency range. This substantiates the assumption made in 
the previous section to reach the Fourier heat law in a generalized sense. 
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IV. Wavevector Dependence of Specific Heat 

In the previous section, we studied the specific heat and other quantities like longitudinal 
viscosity and modulus using a schematic model to show the self-consistency of the relation 
between structural relaxation and the frequency dependence of specific heat. Here we con- 
sider the wavevector and frequency dependent specific heat in a liquid. Starting from the 
generalized hydrodynamic equations for the conserved densities in q-space, we obtain an 
equation, 

iujp c v 5T(q, uj) = -q 2 n5T(q, uj) + - r ( |^ j 5T(q, uj) (15) 

(p uj 2 - q 2 K T (q,uj) \oT J p 

which describes dynamics of the energy fluctuations over different length and time scales. 

Kt(q,uj) is the wave vector and frequency dependent bulk modulus given by, 

K T (q) 

where T(q,co) is the wavevector and frequency dependent longitudinal viscosity devided by 
the square of the speed of sound c 2 (q) = K T (q)/p Q . K T (q) is the zero frequency limit 
of KT(q,uj). The energy equation (|TJ) reduces to the wavevector dependent Fourier heat 
equation, 

tuj5T(q, u) = -q 2 x(q, v)$T(q, uj) (17) 

where x(Q7 u ) = K l i.PoC P (q,uj)) is thermal diffusivity and c p (q, uj) is the g-dependent specific 
heat given by, 

Cp(q,u) = 



1 + (7, - 1) ' 



1 + iujT(q, uj) 

and 7 g is the ratio c p (q)/c v . Here in obtaining the Fourier heat equation (|TTD, we have 
assumed that the quantity, 

M(q, uj)[uj + Li/{q)M{q, uj) 
where u(q) = c 2 (q)/p Q and M(q,uj) = 7(g) + mjT(q,uj). is much smaller as compared to 
unity. In the long wavelength limit this quantity reduces to A(u) given by Eq. (|TT|) . 



V. Results from the Mode coupling Theory 

In this section, we consider the time dependent longitudinal viscosity obtained from a mi- 
croscopic theory of Statistical Mechanics, instead of taking inputs from experimental results 
as was done in the section III. We predict structural aspects i.e. the wave vector dependence 
in the specific heat at different frequencies. 

In the simplest form the self-consistent mode coupling theory predicts a sharp transition 
of the supercooled liquid to nonergodic phase. In later versions it was shown that due to 
coupling of density fluctuations with currents this sharp transition is eliminated - the full 
model with the cutoff mechanism included is termed as the extended model. In the work we 
consider the extended model where the cutoff function is adjusted to obtain agreement with 
the viscosity of the supercooled Liquid to the results obtained from Computer simulations. 
The details of the model and the scheme for computation of the density correlation function 
using the proper cutoff function is presented elsewhere. We consider a one component 
Lennard- Jones system for computing the structural relaxation properties using the MCT. 
The temperature T* and density p* are expressed in the standard units of e/K B and cr 3 
respectively, e is the unit of energy in a L-J system and a is the diameter of a particle. 

In computing the dynamical behavior of the density correlation function we estimate the 
cutoff parameters of the theory so that the shear viscosity obtained from the self consistent 



MCT agrees with computer simulation results 14]. For the simulation results on one compo- 



nent model we use the recent results of Rucco et. al. flTJJ using special techniques that avoid 
the typical problem of crystallization in one component systems. From the self consistent 
results for the density correlation functions we compute the mode coupling integrals for the 
longitudinal part of the memory function related to the decay of the density correlation 
functions. 

The longitudinal viscosity in the zero wave number limit is shown in Fig. 7 for the temper- 
ature range around T c . The longitudinal viscosity shown for the temperature range around 
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T c /T less than 1, follows the power law behavior. The viscosity diverges with exponent equal 
to 1.9 around T c and for lower temperatures ( T c /T ) 1 ) the behavior follows a Vogel Fulcher 
form. This is usual with results of extended MCT |jl6| . We then use this frequency depen- 
dent memory function or the longitudinal viscosity to compute the corresponding frequency 
dependent specific heat. In the Mode-Coupling approximation, the normalized longitudinal 
viscosity, T(q,u), is given by, 



where ip(q, t) and S(q) represent density correlation function and the structure factor respec- 
tively. c(q) is the direct correlation function of the system, q denotes the unit vector along 
q and n is the number density. Using the above expression for the longitudinal viscosity in 



the frequency. This is shown here for three different wave vectors, q<J=0, 7. 05 (peak of the 
structure factor) and 30 (upper cutoff taken for the k-integral) at temperature T*=.559. The 
inset of the corresponding figures shows the secondary peak predicted for fast processes at 
very high frequency window - representing the so called /3-processes in the supercooled liquid. 
The peak in the imaginary part shifts to lower frequency with lowering of temperature. In 
Fig. 9(a) we show the variation of the peak position with temperature of the liquid. The solid 
line in the figure shows Vogel-Fulcher fit with T o =.014. In order to indicate the structural 
dependence we also show in Fig. 9(b), the dependence of the peak position on the wave 
number q at a fixed temperature T*=.559 . The peak frequency signifying the dominant 
time scale for relaxation at different wavenumbers follows the nature of the structure factor. 
It shows a minimum at g-value which corresponds to the peak in the structure factor of the 
liquid. Successive minima in the figure correspond to the other less pronounced maxima in 
the structure factor. 




(20) 



Eq. (P~S|) , we calculate the g-dependent specific heat in the supercooled liquid. 

In figure 8(a) and (b) we show respectively the real and imaginary part of c*(q, u) vs. 
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VI. Discussion 



In section III, we considered two types of experimental measurements on the supercooled 
liquids, respectively related to the energy fluctuations and the structural relaxations. This 
was done to check the consistency of formula for specific heat obtained from simple analysis 
of the Hydrodynamic equations. The comparisons done in section III indicates that the 
frequency dependent specific heat can be understood in terms of the structural relaxation 
data in terms of the analysis proposed in Ref. 0. 

Subsequently we apply the standard forms of the self-consistent Mode Coupling Theory 
to compute the frequency dependence of the viscosity and compute those for the specific 
heat. Here we use the formula in the terms of the Generalized Hydrodynamics, extending 
the model to large k, or small wavelengths. We find that the dispersion in the specific heat 
decreases as we go to the smaller length scales(higher q*) with corresponding increase in 
spectrum width. This demonstrates the fact that at short length scales, the relaxation is 
fast and here the memory effects reflecting cooperativity are not strong. 

The peak position (u p ) in the imaginary part of the c*(q, uj) shifts to higher frequency 
with the increase of the corresponding wavevector q*. It however reaches to a minimum 
frequency at the structure factor peak. Finally since the MCT relates to the two step 
relaxation process in supercooled liquids, there is a corresponding implication on the specific 
heat curve predicting a peak at very high frequency in the specific heat. This is shown in 
the inset of figure 8(b). It is a consequence of secondary relaxation in the supercooled liquid. 
Due to the constraints on the MCT at very low temperatures, we could not study the thermal 
response of the system close to the glass transition temperature T g . However, as is shown 
in the Fig. 9(a), the main peak in the specific heat moves towards the smaller frequencies 
with decreasing temperature, thus at the temperatures very close to the glass transition one 
can expect the two peaks to lie further apart from each other. We have ignored here effects 



of nonlinearities in the energy equation [[L7[ [18fl . This can produce frequency dependence on 
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other transport coefficients like thermal conductivity as well. However, observation of such 
a peak will further strengthen the validity of the simple analysis presented here in energy 
transport in terms of structural relaxation behavior. 
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Figure Captions 

Fig. 1(a): Theory fit to the experimental data(filled circles) for the real part cf p (u) of 
the specific heat (c p (u)) at three temperatures T— 201.4°X(continuous line), 203.9°i ; r(long 
dashed) and 211.4°X(dotted). u>* = ut , where r Q is the units of time used(see text). 

Fig. 1(b): imaginary part c" p (uj) of specific heat corresponding to the real part shown in 
Fig. 1(a). arrows along the frequency axis indicate the peak position in the imaginary part 
of the corresponding viscosity. 

Fig. 2: temperature variation of viscosity in supercooled glycerol as obtained from the 
specific heat fitting. Viscosity is given in units of K r . experimental results for the viscosity 
are shown in the inset. The continuous line( both in the main figure and the inset) is the 
Vogal-Fulcher, 77 = rj exp[A/(T -T c )\, fit with r? o =1.08x 10" 4 P, A= 2480°K and T =128°K. 

Fig. 3(a): real part, M'(u), of the normalized longitudinal modulus M(u;)(see text) is 
plotted at three different temperatures T=203.9°i ; r(solid line), 211.4°i^(dashed line) and 
221.5°i^(dotted line), dots are the experimental results of ref. fpL3| . frequency axis is scaled 
with respect to the peak values u p for the three different temperatures. 
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Fig. 3(b): Imaginary parts of the normalized longitudinal modulus M{uj) corresponding 
to the real parts shown in the figure 3(a). 

Fig. 4: peak position (uj p ) in the imaginary part of the specific heat is plotted as a 
function of the temperature (T). continuous line is the VF fit :uj p = u exp[—A/(T — T c )]; 
with u = 1.0 x 10 15 Hz, £=2559.35°fT and T Q = 128.22°^. 

Fig. 5: Imaginary part of the specific heat, d p {u) at T=21A°K. Arrow along the frequency 
axis at 2.45 indicates the peak position in the imaginary part of the corresponding viscosity 
shown in the inset. 

Fig. 6 real and imaginary parts of A(cj)(see text), for supercooled glycerol, is plotted at 
T = 201.4K. 

Fig. 7: mode-coupling viscosity is plotted as a function of temperature, it shows 
a vogul-Fulcher behaviour, r} exp[A/(T - T Q )] for T < T c with T o =.023 while for higher 
temperatures(T > T c ) it follows a power low behaviour with exponent 1.9 . Arrow along the 
temperature axis at T c /T— 1.36 indicates the power- low divergence. 

Fig. 8(a): MCT results for the real part of the normalized specific heat c*(q,u) = 
(c p (q,uj) — c v )/(c p (q) — c v ) for three wave vectors q*=0 (dotted), 7. 05 (continuous) and 
30(dashed) at T*=.559. Here frequency uo* is in the units of the inverse of Lenard- Jones 
time, r = ^f. 

Fig. 8(b): imaginary part of c p (q,uj) from the MCT corresponding to the real parts 
shown in Fig 8(a). Inset shows the secondary peaks predicted by the MCT for the same 
three q- values. 

Fig. 9(a): variation of the peak position (cu p ) in the specific heat c p (q, uj) with temperature 
for q—0. along the y-axis, we have shown uo* = lot x 10 2 . 

Fig. 9(b): variation of the peak frequency (u p ) with wave vector q*. u p reaches to 
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minimum at g*=7.05 at which the structure factor shows a maximum, uo* = lot x 10 2 . 
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